Vortex Lattice in Planar Bose-Einstein Condensates with Dipolar Interactions 
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In this letter we investigate the effects of dipole-dipole interactions on the vortex lattices in fast 
rotating Bose-Einstein condensates. For single planar condensate, we show that the triangular lattice 
structure will be unfavorable when the s-wave interaction is attractive and exceeds a critical value. 
It will first change to a square lattice, and then become more and more flat with the increase of 
s-wave attraction, until the collapse of the condensate. For an array of coupled planar condensates, 
we discuss how the dipole-dipole interactions between neighboring condensates compete with the 
quantum tunneling processes, which affects the relative displacement of two neighboring vortex 
lattices and leads to the loss of phase coherence between different condensates. 
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Introduction. During the latest ten years an extensive 
study of the many-body physics of ultracold atomic gases 
once again reveals that the interactions between parti- 
cles play a significant role in a many-body system,such 
as determining the structure of a vortex lattice 0, 0, 
and driving quantum phase transitions 0, 0. Recently 
the Bose-Einstein condensate (BEC) of chromium has 
been realized, where the magnetic dipole-dipole inter- 
action is 36 times larger than the alkali atoms[6j; and 
it is very hopeful that the quantum degeneracy of po- 
lar molecules, which have large electric dipole moments, 
may also be achieved soon0- In the study of these sys- 
tems, it is necessary to incorporate the dipolar interac- 
tions between atoms (molecules). Theoretical works have 
revealed many distinctive effects caused by the dipolar 
interactions. For examples, owing to its anisotropic na- 
ture, the stability of a dipolar BEC depends on the trap 
geometry |8|,|9j and its excitation spectrum exhibits roton- 
maxon structure^!; an d owing to its long-range nature, 
the dipolar Bosons in optical lattice provide supersolid 
and checkerboard phases because of the non-negligible 
interactions between nearest-neighboring sites 

In this letter we are interested in how these character- 
istics of the dipolar interactions manifest themselves in 
the structure of vortex lattice in a fast rotating BEC. We 
consider a BEC of chromium which is fast rotated with 
respect to z direction and loaded into a one-dimensional 
lattice imposed along the z direction, as schematically 
illustrated in Fig.JUi). In each site, it is a planar con- 
densate containing a large number of vortices. When 
the sites are widely separated, these condensates are iso- 
lated from each other, and when the sites get closer, they 
are coupled with their neighborhoods via quantum tun- 
nelling and the long-range dipolar interactions. Here we 
consider the case that the spins of atoms are all polarized 
and the spin polarizations are antiferromagnetically ar- 
ranged among the lattice sites, therefore both the intra- 
layer and the inter-layer dipolar interactions are repul- 
sive. This antiferromagnetic arrangement can be realized 
by manipulating the configuration of external magnetic 
field, for example, by usin g th e magnetic field nearby 
permanent-magnet material 



This work contains two parts. First, we consider an iso- 
lated planar condensate under fast rotation. When the 
rotation frequency is very close to the trapping frequency, 
the size of vortex core becomes comparable to the spacing 
between vortices, and its effect becomes important 0, 
Therefore the repulsive s-wave interaction favors a trian- 
gular vortex lattice where the vortex cores are arranged in 
a most symmetric way so that the modulation of the con- 
densate density is spatially homogenous, while the attrac- 
tive interaction favors the lattice structure such as the 
square lattice and the flat lattice, which are less symmet- 
ric than the triangular lattice. However, in the absence 
of the dipolar interactions, such lattices can not exist 
because the condensate itself collapses when the s-wave 
interaction becomes attractive. The presence of dipo- 
lar interaction prevents the condensate from collapse for 
small s-wave attraction and thus the effect of attractive 
interaction on the lattice structure becomes observable. 
Here we are going to show that a triangular lattice will 
be unfavorable in the parameter region that the strength 
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FIG. 1: (color online) (a) A schematic of the model. The 
blank line stands for the one-dimensional optical lattice, and 
the red arrows denote the directions of spin polarization, (b) 
The density profile of a flat vortex lattice, the blank regime 
is the low density regime and the bright regime is the high 
density regime, (c) A schematic of the coincident phase (left) 
and the staggered phase (right). 
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of s-wave attraction exceeds a critical value but is not 
enough to induce collapse. With the increase of s-wave 
attraction, the vortex lattice will experience a structure 
transition from triangular to square, and then to very flat 
lattice as shown in Fig. ^b). 

Secondly, we study the effects of the coupling between 
the neighboring condensates. In the situation that the 
quantum tunneling is dominative, the condensates in dif- 
ferent sites are coherent, and the positions of vortices 
in different layers are coincident with each other (see 
the left one of Fig^c)). However, because the inter- 
layer dipolar interacting energy will be reduced when the 
high density region of one condensate coincides with the 
low density region of its neighboring condensate, it fa- 
vors that the vortex lattices are stagger-arranged among 
the lattices(see the right one of Fig^c)). Here we will 
show that the competition between the quantum tunnel- 
ing and the inter-layer dipolar interactions will lead to 
a quantum phase transition from the coincident phase 
to the staggered phase. Moreover, the tunneling ampli- 
tude will be suppressed exponentially when the vortex 
lattices are staggered, and consequently the fluctuation 
of relative phases between neighboring condensates will 
be enhanced. Therefore, this quantum phase transition 
is characterized by the loss of phase coherence between 
sites, in analogy with the superfluid-Mott insulator tran- 
sition. 

Model and Method. We consider a two-dimensional 
BEC, which is trapped by a harmonic potential mw 2 r 2 /2 
in the xy plane and fast rotated with fixed frequency f2, 
the system is described by following energy functional 

K = J d 2 r^*(r)H y(r) + J d 2 rd 2 r' p(r)V(r - r')p(r'), 

(1) 

where ^(r) is the macroscopic wave function of the con- 
densate, p(r) is the condensate density, and V(r — r') is 
the interactions between atoms which include both the s- 
wave contact interaction and the magnetic dipolar inter- 
action. The single particle Hamiltonian Hq is written as 
— h 2 V 2 / (2m)+mu} 2 r 2 /2— flL, whose eigenstates are Lan- 
dau levels in the fast rotating limit where is very close 
to to, and the interaction energy is much weaker than the 
spacing between different Landau levels 2MI when the 
condensate density is very low, therefore we can restrict 
ourselves in the lowest Landau level(LLL)P|. When the 
vortices form a uniform lattice, the wave function \&(r), 
excluding the Gaussian factor, is a double-periodic ana- 
lytical function, and can be uniquely given by the Jacobi 
theta function|l4|. Although the Gaussian profile of the 
global density is unstable with respect to very weak lat- 
tice distortion^, the LLL wave function serves very well 
as a trial wave function to determine the vortex lattice 
structure. For examples, this wave function has been suc- 
cessfully used to explain the dynamic formation of vor- 
tex stripe observed in a recent experiment [l3l lf l| ; and 



to predict the structure transition of vortex lattice from 
triangular to square in a two-component condensate |l4|. 
which has been observed later 

By denoting bi and X>2 as the basis vectors of the lat- 
tice, a uniform vortex lattice in the LLL is characterized 
by two parameters, which are b2/bi = u + iv describing 
the lattice type and v c = b 2 v denoting the area of a unit 
cell. The condensate density takes the form as p(r) = 

1 / ( 7rcr2 ) Ek 9k exp (iK • r ) exp( -r 2 /cr 2 ) [I El , where the 
summation is taken over all reciprocal lattice vectors 
TOiKi+m2K2, with Ki and K2 being the basis of the re- 
ciprocal lattice. Here a is the condensate radius which is 

given by ^/[aj 2 — irvc 1 ]^ 1 where a± = ^/(mw). The 
coefficients denotes <?k/<?o, an d <?k is explicitly writ- 
ten as g K = (-i) m i+™2+m im2e -^|^V^y^72 ; with 
v c \K\ 2 = (27r) 2 w^ 1 [(umi) 2 + (m 2 + umi) 2 ]. 

In the LLL region the mean value of the single particle 
Hamiltonian Eq turns out to be h(oj^ — £l)cr 2 |l4|], which is 
independent of the vortex lattice structure. The energy 
of s-wave interaction, V s (r — r') = gS(r — r'Lcan also 
be evaluated as E s = g J d 2 rp 2 (v) ~ g//(7rcr 2 )[lj], where 
I denotes the summation of J2k I^k/soI 2 - When the s- 
wave interaction is repulsive, i.e. g > 0, we can find that 
the triangular lattice is stable by minimizing E s , and the 
minimization of the total energy, including Eq and E s , 
yields a finite value of condensate radius a. However, 
when g is negative, the minimization of the total energy 
results in a vanishing, which implies the collapse of the 
condensate. 

Single Planar Condensate. Within the LLL mean field 
ansatz, we are going to investigate how the dipolar inter- 
actions affect the structure of vortex lattice of a planar 
condensate. Considering a very deep lattice, at each lat- 
tice site the motion of atoms along z direction is strongly 
confined, and the magnetic dipole moment p of atoms 
are also polarized along z axis, the magnetic dipolar in- 
teraction in xy plane reads 

V d (r 1 ,r 2 ) = M f 2 1 (2) 
4-7T |ri — r 2 | a 

and the dipolar interaction energy E^ is given by 
f d 2 r 1 c? 2 r 2 p(r 1 )Vd(r 1 — r 2 )p(r 2 ). Denoting the relative 
displacement r = Yi — r 2 , and the center of mass dis- 
placement R = (ri + r2)/2, we can first integrate R out 
and Ed turns out to be 




This integral equals to 2vr J dr J (Kr)e- r2 /(2a2) /r 2 , 
which has shot-range divergence and therefore requires 
a short-range cut-off A. Here A should be choose as 4a z 
where a, is the longitudinal harmonic length of each lat- 
tice sitejl 6|. It is important to notice that a dimension- 
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FIG. 2: (a): The dimensionless value of E, which is defined 
as al + W, as a function of the lattice angle 9 with a — 1.92, 
where 6 = arccos (bi ■ b2)/(bi&2)- The insert is a zoom-in 
plot of the range 8 £ [tt/3,tt/2]. (b): £ as a function of 
r = |t>2 /bi | with a — 1.60. (c) the lattice angle of the ground 
state # m in as a function of parameter a and (d) the radio r of 
the ground state as a function of a. Here C denotes collapse, 
F denotes flat lattice, S denotes square lattice and T denoting 
triangular lattice. 
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with r — r/aj_, and the remainder term depending on 
A is proportional to /, thus the total dipolar energy of 
a rotating condensate can be viewed as an isotropic and 
locally repulsive part plus an attractive part, because W 
is always negative. The total interaction energy including 
the dipolar and s-wave interaction turns out to be 
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where a denotes aj_/ A + 4ga±/ (/io/i 2 ) which can be grad- 
ually changed in a wide range because the s-wave inter- 
action constant q can be tuned via Feshbach resonance 
technique^. The structure of vortex lattice in equilib- 
rium can be obtained by minimizing Eun, and the col- 
lapse of condensate occurs when E- lnt becomes negative. 

The results of the minimization of the mean field en- 
ergy are shown in Fig.©. In the regime that the repul- 
sive interaction is dominant, the ground state is a trian- 
gular lattice with 9 mm = tt/3 and |t>2 /bi = 1- When a 



decreases below 2.0, the triangular lattice becomes un- 
stable and will be replaced by the square lattice, which 
can be seen from Fig.JS^) where the energy at 9 = tt/2 
is lower than 9 = 7r/3. Furthermore, the energy takes 
its minimum at the place where |t»2 /t»i is smaller than 
unity when a is smaller than 1.9, as shown in Fig.JSjs). 
From Fig.(|2:)and J2Ji) one can see that the vortex lat- 
tice changes from a triangular one to a square one, and 
then becomes more and more flat until the collapse of the 
whole condensate, as the s-wave attraction increases. In 
this calculation we take A/a± =0.1 which is a typical 
values of current experiments, for chromium the square 
lattice occurs when a s is between — 2.32a and — 2.28ao, 
and the flat lattice occurs when a s is between — 2.55a 
and — 2.32ao, where do = 0.0529nm. 

Coupled Planar Condensates. Now we are going to 
study the effects of coupling between nest-neighboring 
condensates. Provided that the energy scale of the inter- 
layer coupling processes is much smaller than the intra- 
layer energy scale, the structure of vortex lattice will not 
be affected by the inter-layer coupling. Here we are in- 
terested in how the inter-layer dipolar interactions affect 
the relative displacement of vortex lattices in different 
layers, and therefore we focus on a typical case of trian- 
gular vortex lattice. Neglecting the s-wave interaction 
between layers, the inter-layer interaction is 



V d '(ri,r 2 ) 



2d 2 



4tt (d 2 +r 2 f/ 2 ' 



(6) 



where d is the distance between two layers, and r = ri — 
T2 denotes the relative displacement in the xy plane. In 
the large vortex number limit the inter-layer interaction 
energy turns out to be 
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3 jK - ro F(|K|,d) 



where F(|K|,g?) denotes a dimensionless integral 
F(\K\,d) = I dfJ (\K\rd)e- f2d2 /^ f (2 ~ ^ 



(1 



?2\5/2 



(7) 



(8) 



with f = r/d, and ro is the relative displacement between 
two neighboring vortex lattices. Because |(?k| exponen- 
tially decreases as the increase of |K| 2 , it is sufficient to 
keep the terms with small value of \K\ in the summation, 
namely, K = and K x = |K±i l0 | = |K 0) ±i| = |K±i, T i| 
for the triangular lattice. Thus 



d 4ira 2 d 



F(K ,d) + 2CF(K 1 



, d) cos(K 



ro) -(9) 



Here both F(Ko,d) and F(Ki,d) are positive, the con- 
stant C denotes the value of |<7ki /so| 2 (= |5k i/#o| 2 = 
|<?Ki „ i / <7o 1 2 ) which equals to 0.0532 for the triangular 
lattice, and ^ denotes the summation over Ki.o,Kq.i 



4 




0.0 0.1 0.2 0.3 0.4 0.5 

X 



X 



0.5 




(c) 


0.4 






0.3 






0.2 
0.1 


coincident 


stagger 


0.0 













0.0 0.1 0.2 0.3 0.4 0.5 

x 

0.35 x (d) 



CO. 



0.10 0.15 0.20 0.25 0.30 0.35 




FIG. 3: (a,b) The energy E cou (in the unit of i/2) as a 
function of rn = a;(bi + b2), for two typical parameters of 
13. p = 0.1 and E cou takes its minimum at x = for (a); 
P = 0.5 and E cou takes it minimum at x — 1/3 for (b). (c) 
The value of Xmin, which denotes the value where E cou takes 
its minimum, as a function of /3. (d) /3 as a function of lattice 
spacing d, with t fixed at iia[i/(10-Ka 2 a±). 



and Ki._i. Taking the symmetry of the triangular lat- 
tice into account, we set r = x(bi+b 2 ) with x e [0, 1/2], 
E' d takes its minimum at x = 1/3, which means that the 
vortex cores of one lattice displace to the center of unity 
cells of its neighboring lattice. 

Another inter-layer process is the quantum tunneling 
between neighboring condensates. Within the LLL mean 
field theory, the tunnelling energy can be expressed as 
E t = —tcos(f>ex-p(—N v irr 2 ,/v c )/2 where N v is the num- 
ber of vortices, t is the hopping amplitude between two 
neighboring sites, and (f> is the relative phase between two 
neighboring condensates. This term favors ro = 0, which 
corresponds to binding the vortices in different layers at 
same positions in the xy plane fl8L Il9| . Therefore the 
total energy of the processes coupling two neighboring 
condensates, E con , includes E' d and E t . By minimizing 
E cou , it can be found that the competition between these 
two processes will result in a transition of the ground 
state. Two typical curves of -E^ou are shown in Fig|3fa) 
and (b), in the case (a) the vortex lattices are coincident 
and in the case (b) they are staggered. As illustrated in 
Fig|3[c), this transition is driven by the parameter (3 de- 
fined as [io/j, 2 CF(Ki,d)/(tiT<j 2 d), and it is very a sharp 
one. Experimentally, (3 can be gradually increased by 
decreasing the lattice spacing d, with the barrier height 
simultaneously increased to keep t unchanged. As illus- 
trated in Figlsjd), a transition to staggered phase occurs 
when d/a± is decreased below 0.6. 

It is also worthwhile to consider the relative phase fluc- 
tuations between two neighboring condensates. Because 
8(f) is proportional to 1/yi exp(— N v irr 2 /v c ), the denomi- 



nator of which will be suppressed to a exponentially small 
value when r 2 . becomes comparable to v c in the staggered 
phase, the transition of ro is accompanied by the loss of 
phase coherence between different layers and the suppres- 
sion of particle number fluctuations in each site. In this 
sense, this transition is similar to the superfluid-Mott in- 
sulator transition. However, a peculiar point is that this 
transition is driven by nearest-neighboring interactions 
instead of the on-site interactions. 
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Note added. After finishing an initial version of this 
paper, we became aware of the work by Cooper, Rezayi 
and Simon in which the vortex lattice stucture of a planar 
dipolar condensate was also discussed 20]. 
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